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Abstract The kinematic properties of coronal mass ejections (CMEs) suffer from the 
projection effects, and it is expected that the real velocity should be larger and the real 
angular width should be smaller than the apparent values. Several attempts have been 
tried to correct the projection effects, which however led to a too large average velocity 
probably due to the biased choice of the CME events. In order to estimate the overall 
influence of the projection effects on the kinematic properties of the CMEs, we perform a 
forward modeling of the real distributions of the CME properties, such as the velocity, the 
angular width, and the latitude, by requiring their projected distributions to best match the 
observations. Such a matching is conducted by Monte Carlo simulations. According to the 
derived real distributions, it is found that (1) the average real velocity of all non-full-halo 
CMEs is about 514 km s _1 , and the average real angular width is about 33°, in contrast 
to the corresponding apparent values of 418 km s~ 1 and 42.7° in observations; (2) For 
the CMEs with the angular width in the range of 20° — 120°, the average real velocity is 
510 km s _1 and the average real angular width is 43.4°, in contrast to the corresponding 
apparent values of 392 km s _1 and 52° in observations. 
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1 INTRODUCTION 

Expulsions of dense clouds of plasma from the solar corona, which are termed coronal mass ejections 
(CMEs), have been explored in the past four decades. Most CMEs originate in the low er corona, and 
expel about 10 14 - 10 16 grams of plasma with the kinetic energy about 10 30 to 10 32 ergs dHudson et all 



120061) . and the velocity ranging from tens to more than three thousand km s 1 dSt.Cvr et aL , 2000l) . The 



CME phenomenon has provoked great interests within the communities of solar and solar-terrestrial 
physics since CMEs are often associated with ot her activities , such as flares , filament eruptions, radio 
bursts, and solar energetic particle events ( e.g., iMunro eta! [ 1 19791: iKahlerl 11992b IChen etafl. 120061: 
iGopalswamv et all 120081: IChen et all 120091) . Whether a CME can approach the Earth and induce a 
strong geomagnetic storm depends on various properties of the CME, e.g., the location of the source 
region, the angular width, mass, velocity, magnetic field intensity, and so on. 

Since the fir st detection on 1971 December 14 by the seventh Orbiting Solar Observatory (OSO- 
7) coronagrap h ( Tousev et all 1 19731) . CMEs have been observed by several other space instruments 



dYashiro et ail l2008b . The Large angle and Spectrometric Coronagraph (LASCO) on board the Solar 
and Heliospheric Observatory (SOHO) mission has made a significant impact on the understanding of 
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CMEs since its first light in 1996 January. The total number of CMEs that LASCO observed is over 14 
thousand until the end of 2009 February. The huge database, e.g., the CM E catalog in NASA/C DAWfl 
is very useful for investigating various properties of CMEs. For example, lYashiro et al.l d2004l) found 
that the average velocity of normal CMEs (with the angular widths between 20° — 120°) is about 428 
km s _1 , and the corresponding average angular width is 56° for the CME events observed during 1996- 
2002. Although the statistical results might change more or less as more observations are made or 
different methods are applied in identifying the CME events, research on the CME kinematic properties 
based on this catalog should have a high degree of reliability. However, the measurements of the CME 
kinematic properties directly from coronagraph observations generally suffer from the projection effects. 
For example, it would be expected that t he real CME velocity could be higher and the real angular width 



One straightforward method to get the real kinema tic properties of CMEs is to study the limb CME 
events, which are nearly free of the projection effects. iBurkepile et alj d2004l) identify a sample of 111 
"limb" events among 1351 CMEs observed by the Solar Maximum Mission (SMM) in 1980 and 1984- 
1989. The limb events have an average velocity of 519 km s _1 and an average angular width of 52.3°, in 
contrast to 349 km s _1 and 46° for all CMEs observed by SMM. It is understandable that the velocity of 
the limb events is larger than that of all events, however, the result of the angular width is not consistent 
with the projection effects, which might be du e to the limb sample i s biased to the eruptive filament 
association and therefore is not representative dBurkepile et all 120041) . A similar survey with a larger 
sample is strongly in need. It is noted that, however, uncertainty may arise in identifying the limb events 
in this kind of survey. If EUV dimmings are used to identify the source region of the CMEs, many slow 
events without evident dimmings would be missed. 

Se veral attempts have be en tried to correct the projection effects for individual events. For halo 
events. iMichalek et al.l d2003l) presented a correction method and found that halo CMEs have an average 
real velocity of 1080 km s _1 and average real a ngular width of 120° , in contrast to 900 km s _1 and 
360° before the correction. For non-halo CMEs, iLeblanc et alj d200ll) proposed a correction method, 
where the locus of the associated flare is used to determine the center of the source region of each CME. 
lYeh et al] d2005l) slightly modified such a method and applied it to 619 CME/flare events where the loci 
of the flares are available in the Solar Geophysical Data Report. They found that the average velocity of 
the sampled CMEs is 792 km s^ 1 , and the average angular width is 59°, compared to 549 km s _1 and 
77° before correction. It is noted that their surveys may be biased to strong events (so that the loci of the 
flares were well documented), which results in a remarkably high average velocity. 

One common requirement for the aforementioned inversion methods is that one needs to determine 
the central position of the CME source region, which is not easy to be determined for a large sample. In 
order to correct the projection effects without the knowledge of the central position of the source region 
for each CME event, we propose a forward modeling technique in this paper, and apply it to derive the 
real velocity and angular width. This paper is divided into four parts. The forward modeling method 
is described in Section [2] and the results are shown in Section [3] which is followed by discussions in 
Section[4] 

2 METHOD 

The directly measurable quantities of CMEs are the apparent velocity (v), the apparent angular width 
(AW), and the measurement position angle (MPA). With a large sample of events, we can get their 
distributions statistically, which are fitted with respective formulae. Our forward modeling method can 
be described as follows: Assuming that the real velocity, real angular width, and the real latitude of 
sample CMEs follow the same functions as their apparent counterparts, but with different parameters, 
we project the sample CMEs into the plane of the sky to get the projected distributions of their kinematic 
quantities. The parameters are adjusted by trial and error until the projected distributions of v, AW, and 
MPA best fit the observed distributions. 



http://cdaw.gsfc . nasa.gov/CME_ list 
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In order to get the statistical results of the observations, we appeal for the NASA/CDAW CME 
catalog, which is based on SOHO/LASCO observations. There are 14120 CMEs from 1996 January 1 1 
to 2009 February 28 in this catalog. Among these events, 106 CMEs, whose velocity can not be detected, 
and 396 full halo CMEs are excluded, by which we obtain a sample of 13618 CMEs. The distributions of 
their apparent velocity, apparent angular width, and MPA are showed in Figure [1] which can be fit with 
lognormal, Gaussian, and double Gaussian functions, respectively, as shown in Equations ([Til-©. The 
average apparent velocity of the whole sample is about 418 km s" 1 , and the average apparent angular 
width is about 42.7°. 
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Fig. 1 Distributions of the apparent velocity (left panel), the apparent angular width (mid- 
dle panel), and the MPA (right panel) of the sample CMEs. The velocity is binned in 10 
km s _1 intervals, the angular width is in 10° intervals, and the MPA is in 10° intervals. In the 
right panel, the distance of each triangle from the solar disk center represents the fraction of 
CMEs in this direction, and the circle marks the solar limb. 



Since most CMEs propagate radially with an almost fixed angular width, the so-called cone 



model is often applie d to describe the geom e try of C MEs 
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Michalek et all 120031: iBurkepile etafl l2004t ISchwenn et _ 
20091) . With the assumption of the cone model, the kinematic properties of a CME in 3-dimensional 
space are characterized roughly by 4 variables, i.e., the real velocity (v re ), the real angular width 
(AW re ), the latitude (8) and the longitude (</>) of the source center. It is reasonable to assume that 
CME sources are randomly distributed along the longitude in the heliocentric coordinates. For other 
three variables, we assume that the distributions of v re , AW re , and 8 take the same functions as those 
of the apparent velocity, angular width, and the MPA (Equations [TJO, but with modified parameters, 
which are expressed as follows: 
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where fj,, ai, <J2, and 03 are parameters to be determined. 

In order to quantitatively calculate the projection effects on the CME velocity and angular width, 
we perform Monte Carlo simulations, which are described as follows: One million artificial CMEs 
are created, whose real velocity, angular width, and latitude of the central source region satisfy the 
distributions presented in Equations (0]|6]l, and whose longitude of the central source region is randomly 
distributed in the range of 0°-360°. Each of the artificial CMEs is then projected into the plane of the 
sky, through which the distributions of apparent velocity, apparent angular width, and the MPA can be 
derived. The four parameters o\, (?% and 03) are adjusted until the projected distributions reproduce 
the observational counterparts (Equations QJ3), where the similarity is examined with the Kolmogorov 
test. 

It is noted that, when projecting the artificial CMEs into the plane of the sk y, the relation between 
the apparent velocity and the real valu e depends on the choice of the cone models dSchwenn et all 120051: 
IVrsnak et all 120071: bao et~aill2009l) . We assume that CMEs are symmetric, and propagate out radially 
through the corona, and their l eading edge is a part o f sphere, connecting th e cone tangentially like an 
ice-cream, i.e., the Model C in lSchwenn et al.l d2005l) . Besides, according to lZhang et al.1 (12010b . some 
CMEs, especially those near the Sun-Earth line, would be too faint to be detected by coronagraphs 
due to Thomson scattering, which should be taken into account in our study. For that purpose, we 
calculate the white-light intensity of each CME leading edge at 3Rq in the plane of the sky. Those CMEs 

whose white-light intensity falls below the 3<r of the fluctuation level of the solar wind are considered 

I 1 I — — I 

to be missed by coronagraphs (see Zhang et al., 2010, for details). The expression for calculating the 
Thomson-scattered intensity is given by IB tilings! (119661) . In our work we neglect the limb darkening 
effect, so the expression can be simplified as follows, 

1 ri' 2 

I = TraJ n R — / (2Ccos 3 9 - Acos 5 6), (7) 
81 Jo 

where a represents the Thomson scattering cross section, Jq the photosphere radiative intensity at solar 
surface, no the electron density, Rq the solar radius, 8 the angle between the plane of the sky and the 
direction of CME propagation, A = cos 2 flV9-co S * and C = § - V9 ~g os " 8 - (9 ~ co |" 9)3/2 . 



3 RESULTS 

The derived real distributions of the CME velocity, the angular width, and the latitude are expressed as 
follows: 

jpf x 1 r (llltVe-6.09)^ 

F ^ = n to at- exp Z n w ' (8) 

0.53V27TW 2 x 0.53^ 

1 AW 2 

1 9 2 
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Their profiles are depicted in the left column of Figure|2] As these quantities are projected into the 
plane of the sky, the resulting distributions of the apparent velocity, the angular width, and the MPA are 
shown in the right column of Figure|2]as crosses, where they are compared to the LASCO observations 
{triangles). It is seen that the simulated apparent distributions of the CME velocity, the angular width, 
and the MPA match the observations very well, which implies that the derived distributions in Equations 
(|8ifT0T) might well represent the real situation. 
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Fig. 2 left column: The derived real distributions of the CME velocity (fop), the angular 
width {middle), and the latitude 9 {bottom); Right column: The apparent distributions of the 
CME velocity {top), the angular width {middle), and the MPA {bottom), where the triangles 
correspond to the SOHO/LASCO observations, and the crosses to the projection of the real 
distributions. 



According to the derived real distributions of the CME velocity and the angular width { AW re ) for 
all non-full-halo CMEs, the average CME velocity is 514 kins" 1 , which is ~ 23% larger than the 
apparent value, and the average angular width is 33°, which is ~ 23% smaller than the apparent value. 
In partic ular, for those CMEs whose angular widths are between 20° and 120° (the so-called normal 
CMEs in lYashiro et all [2004). the corrected average CME velocity is 510 kms" 1 , which is ~ 30% 
larger than the apparent value (392 km s" 1 ), and the corrected average angular width is 43.4°, which is 
~ 17% smaller than the apparent value (52°). 

It is noticed that the real angular width distribution is peaked near 0°, however, its apparent angular 
width is peaked around 10°, which is consistent with the observations. Our simulations reveal that the 
significant drop of the number of narrow CMEs (with the angular width less than 10°) is due to that 
some of them are two weak to be detected when they propagate not far from the Sun-Earth direction. 
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4 DISCUSSIONS 



The accumulating observations of CMEs enable the statistics of their properties, such as their angular 
size and their kinematics. However, the measurements of the CME geometry and kinematics suffer 
from the projection effects, and it is expected that the real angular width should be smaller, whereas 
the real velocity should be larger. In this paper, with a forward modeling by Monte Carlo simulations, 
we corrected the distributions of the CME velocity and the angular width for the projection effects, 
as revealed by Equations d&HTOl). It is found that (1) for all non-full-halo CMEs, the real velocity is 
averaged at 514 km s" 1 , and the real angular width is averaged at 33°, in contrast to 418 km s _1 and 
42.7° before correction; (2) for CMEs with the apparent angular width between 20° — 120°, the real 
velocity is averaged at 510 km s _1 , and the real angular width is averaged at 43.4°, in contrast to 392 
km s _1 and 52° before correction. The average real velocity of CMEs derived in this p aper is similar to 
that of the limb events, which are free of the projection effects dBurkepile et al.ll2004l) . 

It should be noted that in our forward modeling, it was assumed for simplicity that the three pre- 
scribed variables (i.e., the CME velocity, the angul ar width, and the latitude of the central source region) 
are completely independent. With a large sample, lYeh et all d2005|) did find that ther e is little correla- 
tion be tween the CME velocity and the angular width. However, Yashiro et al.l (120041) and lVrsnak et al.l 
(2007) claimed that wider CMEs tend to be faster. In order to clarify the discrepancy, we collect the 
13618 non-full-halo CME events observed by SOHO/LASCO during the period 1996 January - 2009 
February and display the scatter plot of the relation between the observed CME angular width and the 
velocity in the left panel of Figure [3] It is found that the two quantities present a very weak positive 
correlation, with the correlation coefficient being only 0.23. We also divide the sample into a subsample 
with CMEs observed near solar maximum and another subsample with CMEs near solar minimum, it 
is also revealed that the CME velocity and the angular width present little correlation. For comparison, 
the right panel of Figure|3]displays the scatter plot of the CME velocities and the angular widths of the 
artificial CMEs in our Monte Carlo simulations. It can be seen that the distributions of the scatter points 
in the two panels of Figure [3] are quite similar. The linear Pearson correlation test indicates that the 
similarity between the two panels is as high as 0.95. It is noted that the sparse data points with v > 1000 
km s~ 1 or AW > 120° are not shown in Figure[3]in order to make the figure clear, though they are 
included in our study. 
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Fig. 3 Scatter plots of the relation between the CME velocity and the angular width in obser- 
vation {left) and simulation (right). The solid lines represent the linearly fitted lines. 

In Figure |4] we compare the correlations between the CME velocity and the MPA in both obser- 
vations and simulations. The left panel displays the scatter plot in the observations, where we can see 
that the correlation coefficient between the CME velocity and the MPA is almost zero. The right panel 
shows the corresponding scatter plot in our Monte Carlo simulations. The linear Pearson correlation test 
indicates that the similarity between the distributions in the two panels is as high as 0.91. Similarly to 
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Figure [3] the sparse data points with v > 1000 km s x or MP A > 200° are not shown in Figure |U 
though they are included in our study. 
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Fig. 4 Scatter plots of the relation between the CME velocity and the MPA in observation 
(left) and simulation (right). The solid lines represent the linearly fitted lines. 

It is noted in passing that the CME angular width and the MPA also show little correla- 
tion. Therefore, with the largest-ever sample that includes almost all the CMEs observed by the 
SOHO/LASCO coronagraph, we found that the CME velocity, the CME angular width, and the lat- 
itude of the central source region are all mutually uncorrected. Such a lack of correlation validates 
our assumptions used in our forward modeling in deriving the CME real velocity and angular width 
distributions. 

Acknowledgements We are grateful to C. Fang and M. D. Ding for discussions and suggestions 
throughout the work. This research is supported by the Chinese foundations 2006CB806302 and NSFC 
(10403003, 10933003, and 10673004). The CME catalog used in this paper is generated and maintained 
at the CDAW Data Center by NASA and The Catholic University of America in cooperation with the 
Naval Research Laboratory. SOHO is a project of international cooperation between ESA and NASA. 

References 

Billings, D. W., 1966, In: A Guide to the Solar Corona, New York: Academic Press 
Burkepile, J. T., Hundhausen, A. J., Stanger, A. L., et al., 2004, J. Geophys. Res., 109, A03103 
Chen, A. Q., Chen, P. E, Fang, C, 2006, A&A, 456, 1 153 

Chen, A. Q., & Zong, W. G., 2009, Research in Astron. Astrophys. (RAA), 9, 470 
Fisher, R. R., & Munro, R. H„ 1984, ApJ, 280, 428 
Gao, P. X. & Li, K. I., 2009, MNRAS, 394, 1031 

Gopalswamy, N., Yashiro, S., Michalek, G, et al., 2008, Earth, Moon, and Planets, 104, 295 
Howard, R. A., Michels, D. J., Sheeley, N. R., Jr., Koomen, M. J., 1982, ApJ, 263, L101 
Hudson, H. S., Bougeret, J. -L., Burkepile, J., 2006, Space Sci. Rev., 123, 13 
Kahler, S. W., 1992, ARA&A, 30, 1 13 

Leblanc, Y, Dulk, G. A., Vourlidas, A., Bougeret, J. -L., 2001, J. Geophys. Res., 106, 25, 301 
Michalek, G., Gopalswamy, N., Yashiro, S., 2003, ApJ, 584, 472 
Munro, R. H, et al., 1979, Sol. Phys., 61, 201 

Schwenn, R., Dal Lago, A., Huttunen, E., Gonzalez, W. D., 2005, Ann. Geophys., 23, 1033 

St.Cyr, O. C, et al., 2000, J. Geophys. Res., 105, A8, 18, 169 

Tousey, R., et al., 1973, Sol. Phys., 33, 265 

Vrsnak, B., Sudar, D., Ruzdjak, D., Zic, T, 2007, A&A, 469, 339 

Yashiro, S., Gopalswamy, N., Michalek, G., et al., 2004, J. Geophys. Res., 109, A07105 



8 



Wu, Y. & Chen, P. F. 



Yashiro, S., Michalek, G, Gopalswamy, N„ 2008, Ann. Geophys., 26, 3103 
Yeh, C. T., Ding, M. D., Chen, P. R, 2005, Sol. Phys., 229, 313 

Zhang, Q. M., Guo, Y, Chen, P. R, Ding, M. D., Fang, C, 2010, Research in Astron. Astrophys. (RAA), 10, 461 



